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Abstract 

We study the stabilities and classical solutions of Euler-Poisson equations 
of describing the evolution of the gaseous star in astrophysics. In fact, we 
extend the study the stabilities of Euler-Poisson equations with or without 
frictional damping term to some special dimensional spaces. Besides, by using 
the second inertia function in 2 dimension of Euler-Poisson equations, we 
prove the non-global existence of classical solutions with 2 f^(p \ u\ 2 + 2P)dx < 
gM 2 — e, for any 7. 



1 Introduction 

The evolution of a self-gravitating fluid such as gaseous stars can be formulated by 
the Euler-Poisson equations of the following form: 

p t +V ■ {pu) =0, 
(pu) t +V ■ (pu ®u)+ VP + (3pu =-pV$, 

S t + u ■ VS = 0, U 
A$(t,x) = a(N)gp, 

where a(N) is a constant related to the unit ball in R N : a(l) = 2; a(2) = 2tt; For 
simplicity, we take the constant term g — 1. For N > 3, 

a(N) = N(N - 2)V(N) = N(N - 2) J^ 2 + 

where V(N) is the volume of the unit ball in R N and T is a Gamma function. As 
usual, p = p(t,x), u = u(t,x) G H N and S(t,x) are the density, the velocity and 



*E-mail address: nevetsyuen@hotmail.com 



1 



2 



M.W.Yuen 



the entropy respectively. P = P{p) is the pressure. In the above system, the self- 
gravitational potential field $ = $(f, x) is determined by the density p through the 
Poisson equation. The equations ([T])i and ([I]) 2 are the compressible Euler equation 
with forcing term. The equation (JTJ)3 is the Poisson equation through which the 
gravitational potential is determined by the density distribution of the gas itself. 
Thus, we called the system ([1]) the Euler-Poisson equations. Here, the viscosity term 
does not appear, that is, the viscous effect is neglected. In this case, the equations 
can be viewed as a prefect gas model. For A/" = 3, (jTJ) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [2J, [S] for a detail about the system. 

If we take S(t, x) = hiK, for some fixed K > 0, we have a 7-law on the pressure 
P(p), i.e. 

P{p) = Kp^ = P -, (2) 

7 

which is a commonly the hypothesis. The constant 7 = Cp/cy > 1, where cp, c v 
are the specific heats per unit mass under constant pressure and constant volume 
respectively, is the ratio of the specific heats, that is, the adiabatic exponent in ([2]). 
In particular, the fluid is called isothermal if 7 = 1. 

For the physical dimension iV = 3, we are interested in the hydrostatic equilib- 
rium specified by u = 0, S = InK. According to [2], the ratio between the core 
density p(0) and the mean density p for 6/5 < 7 < 2 is given by 

P [-?> ■ i X 



z=z 



P (o) v z 

where y is the solution of the Lane-Emden equation with n = 1/(7 — 1), 
y + -y + y n = 0, y(Q) = a > 0, y{0) = 0, n ' 



z 7 — 1 

and z is the first zero of y(z ) = 0. We can solve the Lane-Emden equation 
analytically for 



yanal \%) 



1 - \z\ n = 0; 

n = l; 
n = 5, 



i 



and for the other values, only numerical values can be obtained. It can be shown 
that for n < 5, the radius of polytropic models is finite; for n > 5, the radius is 
infinite. 

Gambin [5] and Bezard pQ obtained the existence results about the explicitly 
stationary solution (u = 0) for 7 = 6/5: 

V*^. (3) 
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The Poisson equation ([I]) 3 can be solved as 

$(t,x)= G(x -y)p(t,y)dy, 
Jr n 

where G is the Green's function for the Poisson equation in the iV-dimensional spaces 
defined by 

(\x\, N=l; 
G(x) = < logjac|, iV = 2; 

[^=7, N>3. 

In the following, we always seek solutions in spherical symmetry. Thus, the Poisson 
equation ([I]) 3 is transformed to 

r N - l $ rr (t, x) + (N — 1) r N ~H r =a (N) pr N ~\ 

®r=^ f p(t,s)s N ^ds. 
r Jo 

In this paper, we concern about existence of global solutions for the N- dimensional 
Euler-Poisson equations, which may describe the phenomenon called the core col- 
lapsing in physics. And our aim is to construct a family of such blowup solutions to 
it. Historically in astrophysics, Goldreich and Weber [6] constructed the analytical 
blowup solution (collapsing) solution of the 3-dimensional Euler-Poisson equation 
for 7 = 4/3 for the non- rotating gas spheres. After that, Makino [8] obtained the 
rigorously mathematical proof of the existence. And in [I], we find the extension of 
the above blowup solutions to the case N > 3 and 7 = (2N — 2)/N. In [11], the 
solutions with a form is written as 

r2/( _r_ )W -2) ? forr<a(t)Zp 



p(t,r) = ^ Ny ^> ' r ' " , u(t,r) =$ r , S(t,r) = LnK 

I 0, tor aytjZn < r. 

a ( t ) = ~~ o(t)^-i » a (°) = a o > °' a(0)= °i> 
y{z)+^y{z) + N %-^ y{z) N,{N - 2) = P, y(0) = a > 0, y(0) = 0, 

(4) 

where p = [N(N - 2)X]/(2N - 2)K and the finite is the first zero of y(z). 

Recently, in Yuen's results [12] there are existing a family of the blowup solution 
for the Euler-Poisson equations in radial symmetry in the 2-dimensional case, 

p t + Up r + pU r + \pU = 0, ^ 
p (U t + UU r ) + Kp r = -2^ p(t, S)sds. 

Those are: 



p(t,r) 



a(t) 



1 p y(r/a(t)) 



u(t, r) 



-ait) 1 

a W = _ ' a (°) = a > °> a(0)= ai; 
j^+i^+f e y(x) = p, y(0) = a, y(0) = 0, 



(6) 
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where K > 0, \i — 2X/K with a sufficiently small A and a are constants. 

(1) When A > 0, the solutions blow up in a finite time T; 

(2) When A = 0, if a\ < 0, the solutions blow up at t = —a /ai. 

In this paper, we study the stability of the Euler-Poisson equations with or 
without frictional damping in Section 2. After that, we use the second inertia 
function in 2 dimension, that is 

H{t) = I p(t,x) \x\ 2 dx, 
Jq 

to prove the following result: 

Theorem 1 Let (p(t),u(t)) be a classical solution with the property 

2 [p \u\ 2 + 2P)dx < gM 2 - e, 
Jq 

of the Euler-Poisson equations (TJP in 2 dimension on <t < T , here 

M = p(t,x)dx, 
Jq 

is the total mass which is constant for any classical solution and e is an positive 
uniform constant. Then T must be finite. 

Initially, in [10], Sideris applied the second inertia function to obtain the insta- 
bility result for the Euler equations. The similar result for the Euler equations with 
frictional damping can be found in [7]. The corresponding instability analysis for 
the 3-dimensional spherically symmetric Euler-Poisson Equations for 7 > 4/3 was 
appeared in [9]. Besides in [3], the instability similar result of the 3-dimensional 
Euler-Poisson equations for any bounded domain was obtained later. In a similar 
way, we extend the results for the case when the some special dimensions N > 2. 



2 Stability 

In this section, we study the stability of Euler-Poisson equations in higher dimensions 
(N > 3). The total energy can be defined by, 

{ E(t) = J a (|p \u\ 2 + ^jP+ ip$) dx, for 7 > 1; 
\ E(t) = J n (\p \u\ 2 + Kp\np + |p$) dx, for 7 = 1. 



In [3], the authors calculated the stationary energy for 3-dimensional the Euler- 
Poisson equations. Here, we generate the result for (iV > 2)-dimensional Euler- 
Poisson equations. 
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Lemma 2 Suppose the solutions p have compact support. For 7 > 1 in Euler- 
Poisson equations, if u = 0, the total energy can be expressed by the integral of the 
pressure P. 

(1) For N = 2, we have, 

gM 2 If f 
E{t) = 4 ( 7 _ ^ + 2 J J P(v) lri \x -y\dydx. 

(2) For N > 3, we have, 

2N - 2 - f 

£ W = (,-!)(* -2) 

Proof. Set u — 0, p(x,t) = p(x), $(x,t) = $(x). 
(l)For N = 2, when u = from <^ 2 , 

VP(x) = -p(x)V$(x). (8) 

We have, by the divergence theorem 

2 / P(x)dx = - I x ■ VPdx (9) 
J Jn 

= J x-p{x)V$dx 

Jn Jn \x — y\ 
= gh, 

where, 

h = I P(x) I ^ - V) • j X - V) dydx + / p(x) / ^"^ yds 
Jn Jn \x - y\ Jn Jn \x — y\ 

= M 2 — h. 
Therefore, we have, 

1 9 

2 ' 



and 



The total energy is 



P(x)dx = -M 2 . 
4 



E = I [ \p\u\ 2 + — —j p ^j dx + \ I P® dx 



n 



gM 2 1 
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(2)For N > 3, first we claim the identity, 

N [ Kp{xfdx = -( N ~ 2 ^ I p ( x )§( x )dx. 
Jn 2 J Q 

To prove the claim, 

Jn \x — y\ Jn \x — y\ 

Multiplying this by x and integrating over Q, we get 

/ x ■ p(x)V$(x)dx = g(N - 2) f p(x) [ - V } ' X dydx. (10) 

Jn J Jn \x — y\ 

With the help of divergence theorem, we get 

N [ Kp r dx = - [ x ■ V (ifp 7 ) dx, 
Jn Jn 

and from (jSJ), 



N / Kp r dx = - I x ■ VPdx = / x ■ pVQdx. 
'n Jn Jn 



From ffTUl). we have 



N [ Kp(x) r dx = g(N -2) [ p(x) [ P{v) } X V }' % dydx = g ■ h, (11) 
Jn J Jn \x — y\ 

where, 



A = (N - 2) / „ W / emEzpE dydx 

x-y\ 

(n-2) ( / P (,) / gM(jp)_|^) d!/fa+ / ^-fy iyix 

n \x — y\ J Jn \x — y\ 



< A '- 2 ) M*) / , P{y ] N _ 2 dydx-h. 
'n \x - y\ 



Thus, we have 

L = — z — / o(x) 

'n \x — y\ 



h = N n I p( x ) I ; — P ^ N - 2 d y dx - 



Hence fJTTjl becomes 



2 Jn Jn \x — y\ 2 



N I P(x)dx = — '- I p(x) I H ' ' N _ 2 dydx 



^ p( x ) I , gp v lN -2 d y dx 

\x-y\ 



2 

N-2 
2~ 



n Jn 



p(x)$(x)dx. 



n 
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The total energy of stationary solutions to Euler-Poisson equations is: 

E — [ — - — P(x)dx + - [ p(x)§(x)dx 
Jul- 1 2 Jn 

1 N f 

-P{x)dx — — — - / P{x)dx 



n 7-1 N-2 JQ 
2N - 2 - iV 7 



(7"l)(iV-2) 



P(x)dx. 



n 



The proof is completed. ■ 

Besides, for the Euler-Poisson equations with frictional damping term, we have 
the energy decay estimate lemma: 

Lemma 3 For the Euler-Poisson equations, suppose the solutions (p, u) have com- 
pact support in Q. With N > 3 and [3 > 0, we have, 

E(t) = ~(3 p\u\ 2 dx < 0. 
Jn 

Furthermore, if f3 = 0, then we have, 

E(t) = 0, i.e. E(t) is conserved. 

Proof. Observe 

Q p ' M ' 2 ) = ( pu ^- u -\pt\ u \ 2 ■ ( i2 ) 

Multiplying -\\u\ to (HJ i gives, 

1 2 1 2 

And multiplying u to yields, 

(pu) t ■ tt + it ■ [V ■ (pu®u)]+u- VP = -pV$ -u- [3p \u\ 2 . (14) 
We sum up f|T3|) and (fT4l) to obtain, 





\ 1 1 1 




M 


/ t 2 1 1 



Notice that 

AT N N 

U ■ [V • (/CU (g) «)] = [<9j (pMj) Mi + pUjdjUi] = \u\ 2 <9j (pUi) + pUiUjdjUi. 

i, j=l i=l i, j=l 
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Hence, we have, 



r i 2 

/ — 77 \u\ V ■ (pu) +u - [V • (pu <S> u)\ dx 
in 2 

r \i 2 N N 
' N N 

u ■ diupui + piiiUjdji 

JV V 

pUiUjdiUj + pUiUjd. 



dx 



dx 



dx 



0. 



Besides, the time derivative for the integrand of the second term in (JTJ)i is: 

Pi = K 1P ^ x d t p 

3 



K'yp'^dipui - ^2 KlP 1 ^ 

i=l i=l 
3 3 



U; 



i=l i=l 

Integrating over Q and applying integration by part to it, we have 

3 



/ P t Gfe = (7 — 1) / d-iPuidx. 
Jn i=1 Jn 



Thus, we get, 



1 



P+cfx — I u ■ V Pdx. 



7 - 1 Jn Jn 
Applying the divergence theorem and using (JTJ) lj3 , we have 



pV$ • udx 



52 "(JVjfl 



(15) 



a(JV)<? 



- / p$ t dx. 
In 
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Therefore, from (TT5]) and (JTJ) i , we get 

- / pV$ • udrr = — ~ I ($d t p + pd t $) dx = -- f d t (p$)dx. 
Jn 2 J n 2 J Q 

By taking the integration for the energy (J7|), we have 

E(t) = [ \-p\u\ 2 + —^—P + -p§ ] dx = -/3 [ pu 2 dx<0. 
AA 2 7-1 2 J t Jn 

For 7 = 1, we have, 



and 



p t dx = — J V (pit) dx = J pu ■ ndx = 0, 



— / Kp In pdx 



f^Pt hi pdx + / Kp t dx 
Jn 

K / V (pit) In pdx 



Similar calculation gives 



uVPdx. 



E(t) = -p I p\u\ 2 dx < 0. 



The proof is completed. ■ 

In here, we calculate the second derivative of the second inertia function for the 
Euler-Poisson equations (N > 2). 

Lemma 4 Consider the Euler-Poisson equations without frictional damping, 

p t +V ■ (pu) =0, 
(pu) t +V ■ (pu ®m) + VP =-pV$, 
A$(t,x) =a(N)gp, 
P =Kp\ 

with (p, u) is the solution with compact support in Q, consider the second inertia, 
i.e. 

H(t) = / p(t,x) \x\ 2 dx. (16) 
Jn 

We have 

' H(t) = 2 /„ [(p \u\ 2 + NP)dx + ^p$] dx, for N > 3; 
F(t) = 2/ n (p|u| 2 + 2P)da;-^M 2 , /or N = 2. 
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Proof. Taking derivatives to (|T6l) implies: 



H(t) = — I V • (pu) \x\ dx = 2xpudx, 
In Jn 



H(t) = 2 x [-V • (pu (8 it) — VP - pV$] dx. 



(1* 



We split (1181) into three parts in the following. For the first part, we apply integration 
by part, 



» N N 

/ x ■ [V • (pw <g> w)] dx = — / Xj 9, (pUiUj) dx = 
Jn j=1 Jn i=l Jn 



pu 2 dx. 



For the second part, we have 



x ■ VPdx = / NPdx. 
v. Jn 



If AT > 3, from (TTO]) and (fTI]) . we get 



N-2 . 

— I px ■ \7&dx = — - — / Qpdx. 



Finally, we have 



H(t) = 2 



N — 2 

(pu 2 + NP)dx H — p$ 



dx. 



If = 2, since from (Q, we have, 



• VQdx = — gM 2 

2 



it gives 



#(£) = 2 (p \u\ z + 2P)dx - gM 2 
'n 



The proof is completed. ■ 

By applying the above lemma, we have the following theorem. 



Theorem 5 Suppose (p, u) is a global classical solution in the Euler-Poisson equa- 
tions without frictional damping. We have 
(l)ForN = 3,4, 



lim inf 



lim inf -^7^ > 

R(t)\n\ 



M 



1/2 



( 7 -l)/2 



> 



{N~t-2(N-1))K M ^-x 
7-1 



1/2 



ifl > and E > 0; 

*/7 > ^TT 11 and E = °- 
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with R(t) = max xe sj(i) {\%\}- Here 



M = I p(t,x)dx, 
In 



is the total mass which is constant for any classical solution. 
(2)ForN > A, 

(2a)the solutions blow up in finite time T, if 1 < 7 < 2(iV — 1)/JV and E < 
(2b)the solutions blow up in finite time T , if 1 < 7 < 2(JV — 1)/N and E = 



Proof. (l)For JV = 3, 4, 7 > 2(JV - 1)/JV and E > 0, we get from Lemma ED, 



— — - / p&dx 
v. 2 J Q 

4 — JV 2 JV7 - 2(JV - 1 



H{t) = 2 / [p \u\ 2 + JVP] dx + 



-pM + 



'n L * 

> 2(JV-2)P > 0. 



7-1 



-P 



cJx + (JV - 2)P 



That is 

P(t) > H(0)+H(0)t + (JV - 2)Et A 
On the other hand, we have, 



P(0)+P(0)t + (JV - 2)Et 2 < H (t) < R(t) 2 M. 



That is, 



lim inf — — > 



t— >oo 



(N-2)E 
M 



1/2 



For JV = 3, 4, 7 > 2(iV - 1) /JV and E = 0, we have 

1/7 



M = I pdx<[ I p^dx 
'n \./n 



111 



(7-l)/ 7 



and 



#(t) = 2 



4 -AT JV T -2(JV-1) D 
p|w| H ; P 



7-1 



cJx + (JV - 2)E 



> 2 / ^-^-^Pd,. 

7 -l 



Thus, we get, 



2[JV7-2(JV-1M 

7-1 ' ' 
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P(0)+#(0)t4 



[Nj -2(N-l)]K ^.^^^ TT ^ ^ 



7-1 

1. [Nj — 2(N — 1)] K 



n\ L - J AfV< if(t) < P(i) z M, 
i?(t) 2 M 



7-1 



in 



1-7 



M 7 <- 



This gives 



lim inf 

t— *oo 



R(t) \tt 



(7-l)/2 



> 



(JV 7 - 2(JV - l))-fT 
7^ 



M 7_1 



nl/2 



(2a)For iV > 4, 7 < 2(N - l)/N and £ < 0, we have 

H{t) = 2 j(iV-2)P + 
< 2{N-2)E. 



2 7 — 1 



That is, 



H(t) < H(0) + H(0)t + {N- 2)Et 2 . 



As £ < 0, if (0) + if (0)t + (N- 2)Et 2 will become negative if t is sufficiently large. 
But, H(t) is always non-negative. Thus we have a contradiction. When iV > 4, 
7 < 2(iV — 1)/N and i£ < the classical solutions blow up after some finite time T. 
(2b)For N > 4, 7 < 2(iV - 1)/N and E = 0, we obtain, 



if (0 = 2 <j (N - 2)P + 

iV 7 - 2(N - 1 



4-JV |2 , jV 7 -2(JV-l) „ 
p|w| H ; 1 



< 2 

Form (fl9|) . we get, 

#(*)< 
After solving it, we have, 



2 

Pete 



7 



7-1 

2(Ny - 2(N - 1))K 



n 



1-7 



M 7 . 



7 



i?(t) < tf(0) + F(0)t + 



(* T - 2(^-1))* 
7-1 



As 7 < 2(iV — 1)/N, we use the above same argument in (3a) to prove non-global 
existence of the solution exists. Thus, when N > 4, 7 < 2(N — 1)/N and E = the 
classical solutions blow up after some finite time T. The proof is completed. ■ 

By applying the similar method in the above section to the Euler-Poisson equa- 
tions with frictional damping (3 > 0, 



p t +V • (pu) =0, 
(pu) t +V ■ (pu ® u) + VP+ =-pV$ - /3pu, 
A$(t,x) =a(N)gp, 
P =fTp 7 . 



(20) 
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Theorem 6 Suppose (p, u) is a global solution of the Euler-Poisson equations with 
frictional damping, i. e. lfW\) . for N > 4, 

(1) the solutions blow up in finite time T , if 1 < 7 < 2(iV — 1)/N and E(0) < 0; 

(2) the solutions blow up in finite time T , if 1 < 7 < 2(iV — 1)/N and E(0) < 0. 



Proof. (l)For iV > 4, 7 < 2(N - 1)/JV and E(0) < 0, we get 

A^7 — 2(N — 1 



H(t) 



l H{t) + I 



A — N |2 

o p \u\ 



P 



7 



dx + (N - 2)E 



1 ■ 2(N-2) 



From Lemma [31 we have, 



1 



Hit) < --Hit) + 



2(JV-2)E(0) 



That is 



H {t) <C 1 + C^ + 2iN -f i0 \. 

P 



As E(0) < 0, C 1 + C 2 e- I3t + [2(N - 2)E(Q)t] /(3 will become negative if t is sufficiently 
large. But, H(t) is always non-negative. Thus we have a contradiction. When 
N > 4, 7 < 2(N-1)/N and i?(t) < for t > 0, the classical solutions blow up after 
some finite time T. 

(2)For iV > 4, 7 < 2(N - 1)/N and E(0) = 0, we get 



ir(t) = + ! 



4-iV . ,2 , jV 7 -2(iV-l) D 

Pp| H : P 



7-1 



cfe + (N - 2)E 



AT 7 _ 2{N - 1) 



Pdx. 



7 



This implies 



Since 



if (0 < Ci + C 2 e 



-/ft 



iV 7 _ 2(JV - 1) 
7^ 



Pete 



M 



pete < 



1/7 



111 



we get, 



if (t) ^ Ci + C^' 3 ' + 



p^dx 

2{N^-2{N - l))K 



(7-l)/ 7 



M 7 Ifi 



1-7 



t. 



^(7-1) 

As 7 < 2(N — 1)/N, we use the same argument in (1) to give the non- global existence 
of the classical solution. Thus, when iV > 4, 7 < 2(N-l)/N and J5(t) < for t > 0, 
the classical solutions blow up after some finite time T. The proof is completed. ■ 
By the similar argument, when 7 = 2( - 7 ^7 1 ^ for iV > 4, any perturbation with 
negative total energy leads to non-global existence. 
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3 Non-global Existence for 2-d Classical Solutions 

We present the proof of the theorem [1] here. 

Proof. From the second inertia function in 2 dimension for the Euler-Poisson equa- 
tions (|17|) 9 without frictional damping ((3 = 0), we have 



as H(t) cannot be negative. Thus the classical solutions of such kinds in global 
existence are not possible. The proof is completed. ■ 
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